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The breaking of chiral symmetry of light quarks at zero temperature in presence of strong 
quantizing magnetic field is studied using Nambu-Jona-Lasinio (NJL) model with Thomas-Fermi 
type semi-classical formalism. It is found that the dynamically generated light quark mass can 
never become zero if the Landau levels are populated and the mass of light quarks increases with 
the increase of magnetic field strength. 



1. INTRODUCTION 

o : 

, The theoretical investigation of properties of compact stellar objects in presence of strong quantizing magnetic field 
■ have gotten a new life after the recent discovery of a few magnetars [1~4]. These stellar objects are believed to be 
1— H , the strongly magnetized young neutron stars. The surface magnetic fields are observed to be > lO-'^^G. Then it is 
I ^ ' quite possible that the fields at the core region may go up to lO^^G. The exact source of this strong magnetic field is 
of course yet to be known. These objects are also supposed to be the possible sources of anomalous X-ray and soft 
gamma emissions (AXP and SGR). If the magnetic field is really so strong, in particular at the core region, they must 
affect most of the important physical properties of such stellar objects and also some of the physical processes, e.g., 
^ . the rates / cross-sections of elementary processes, in particular the weak and the electromagnetic decays / reactions 

' taking place at the core region. 
CnI ■ The strong magnetic field affects the equation of state of dense neutron star matter. As a consequence the gross- 
properties of neutron stars [5-8], e.g., mass-radius relation, moment of inertia, rotational frequency etc. should 
change significantly. In the case of compact neutron stars, the phase transition from neutron matter to quark matter 
which may occur at the core region is also affected by strong quantizing magnetic field. It has been shown that a first 
order phase transition initiated by the nucleation of quark matter droplets is absolutely forbidden if the magnetic field 
strength ~ lO^^G at the core region [9,10]. However, a second order phase transition is allowed, provided the magnetic 
field strength < 10^°G. This is of course too high to achieve at the core region. The study of time evolution of nascent 
quark matter, produced at the core region through some higher order phase transition, shows that in presence of 
Q ' strong magnetic field it is absolutely impossible to achieve chemical equilibrium (/3-equilibrium) configuration among 
|H ', the constituents of the quark phase if the magnetic field strength is as low as B„i ^ lO^^G. 

The elementary processes, in particular, the weak and the electromagnetic decays/reactions taking place at the 
core region of a neutron star are strongly affected by such ultra-strong magnetic fields [11,12]. Since the cooling of 
neutron stars are mainly controlled by neutrino/anti-neutrino emission, the presence of strong quantizing magnetic 
k> ■ field should affect the thermal history of strongly magnetized neutron stars. Further, the electrical conductivity of 
, neutron star matter which directly controls the evolution of neutron star magnetic field will also change significantly 

d, [12]. 

Similar to the study of quark- hadron deconfinement transition inside neutron star core in presence of strong quantiz- 
ing magnetic field, a thorough investigations have also been done on the effect of ultra-strong magnetic field on chiral 
symmetry breaking. In those studies, quantum field theoretic formalisms were mainly used [13-19]. In reference [20] 
Inagaki et al have studied the chiral symmetry violation with NJL model using quantum field theoretic approach in 
presence of strong quantizing magnetic field. In many of these papers, the effect of curvature with or without external 
magnetic field on chiral symmetry violation have been investigated. In the studies by Gusynin ct al in [13,21], have 
thoroughly investigated the chiral symmetry breaking in presence of strong external quantizing magnetic field. They 
have used NJL model in 2 -I- 1 and also in 3 -I- 1 dimensions. It has been shown that the external magnetic field acts as 
a catalyst to generate fermion mass dynamically. In the first paper [13] they have studied it in 2 + 1 dimension and 
showed how the external magnetic field generated dynamical mass of fermion and broke the dynamical flavor sym- 
metry. They have further shown by using NJL model that chiral symmetry breaks dynamically even if the attractive 
interaction between the fermions is extremely weak. In the second paper [21] they have extended the calculation to 
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3 + 1 dimension. In the very nice piece of work by [17] Lee et al have stiidied the breaking of chiral symmetry for 
fermions in presence of external magnetic field. It has also been shown in this work that the symmetry is broken 
dynamically and further the effect of finite density and the temperature of the system on the chiral properties of the 
fermions have been investigated thoroughly in this paper in presence of strong magnetic field. It has been reported 
in this paper that there exists a critical density (or chemical potential) above which the chiral symmetry is again 
restored (which actually indicates the restoration of chiral symmetry at high enongh density) and if it is treated as 
a chiral phase transition, the order will be of first order in nature. On the other hand the chiral symmetry is again 
restored at high temperature above some critical value. In this case the transition is of second order in nature. In 
an extensive review work [22], Klevansky has reported the dynamical chiral symmetry breaking in presence of strong 
external quantizing magnetic field using NJL model in SU{2) and SU{3) flavor space. In this paper the eflFect of 
density (i.e., finite chemical potential) and temperature of the system on chiral symmetry restoration have also been 
reviewed. 

In the present chapter we shall study the effect of strong quantizing magnetic field on the chiral properties of 
QCD vacuum with the help of NJL model following a semi-classical Thomas-Fermi type mean field approach in 
presence of strong quantizing external QED magnetic field. Now in NJL model, there is no in-built mechanism of 
color confinement, however, it can produce two chirally distinct phases- appropriate for confined quark matter within 
the bag and the matter outside the bag. These phases are also known as the Wigner phase and spontaneously broken 
chiral phase respectively. Therefore, if one re-formulates the NJL model in presence of strong quantizing magnetic 
field, it is quite possible to obtain the eff'ect of quantizing magnetic field on these two chirally distinct phases and hence 
obtain the effect of magnetic field on chiral symmetry breaking. Further, it is also possible to obtain bag pressure 
from the diflFerence of vacuum energy densities of these two phases and hence its variation with strong magnetic field. 
Assuming that the confinement and spontaneously broken chiral symmetry are synonymous, Bhaduri et. al. obtained 
some estimate of bag constant from the difference of energy densities [23] for the conventional case. In the present 
chapter we shall modify these original calculations of Bhaduri et. al. [23] and Providencia et. al. [24] to study the 
breaking of chiral symmetry of light quarks in presence of strong magnetic fields and show that the chiral symmetry 
always remains broken in presence of strong quantizing magnetic field if the Landau levels for quarks are populated. 
Our motivation in this work was to study the effect of strong quantizing magnetic field on two chirally distinct phases 
and then obtain the vacuum pressure as a function of strong external magnetic field. Unfortunately, we have noticed 
that Wigner phase does not exist if the Landau levels of quarks are populated and in this formalism, there is no way, 
either by controlling chemical potential i.e., the density of matter (which is meaningless in our investigation since 
we have considered QCD vaciium state) or temperatiirc of the system to restore chiral symmetry. Our study is now 
basically an application of the formalism developed recently to study the equation of state of dense fermionic matter 
of astrophysical interest in presence of strong quantizing magnetic field [25]. 



2. BASIC FORMALISM 



We start with the density matrix p{x,x'), defined by 

p(a;,x')- ^ V(^)V^(mA-|p, ]) (1) 

spin,p 

where tp and tp^ are respectively the negative energy Dirac spinor and the corresponding adjoint, satisfy the equation 

htp = E-ij (2) 
(and similarly for tp^) with the single particle Hamiltonian 

h = -f5^.{p-qfA)+(3m (3) 

with 

2=(S\5). W 

75 and /? are the usual Dirac matrices, A is the ultra-violet cut off in the momentum integral over ( since we are 
considering vacuum, unlike a many body fermionic statistical system we have to put the cut off by hand) and A is the 
electromagnetic field three vector corresponding to the external constant magnetic field of strength along ^-axis. 
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Here the light quark mass m is assumed to be generated dynamically. Now in presence of strong quantizing magnetic 
field along z-direction, the up and down spin negative energy spinor solutions are therefore given by 



il}{x) = 



(T,i) 



(5) 



where 



,(T) 



[2£;_(£;_ -m)]V2 



and 



(i) 
v_ = 



[2E_{E_-m)YI^ 









\ {E_-m)I,-i J 



\ 



(6) 



(7) 



where E- = —{jpl + + 2uqfBmY^^ = —E^, is the single particle energy eigen value, v = 0,1, 2..., are the Landau 
quantum numbers, g/ is the magnitude of the charge carried by /th flavor and 



1/4 



1 



(iy!)l/2 



2-^/2 



exp 



1 



-q/Bm X 



Py 
q/Bm 



{qfBrr^y/^X 



Py 
qfB„ 



(8) 



with H„ is the well known Hermite polynomial of order p, and Ly, are respectively length scales along Y and Z 
directions. Now it can very easily be shown that u = state is singly degenerate, whereas all other states are doubly 
degenerate. We now express the density matrix, as the modified version of Wigner transform in presence of strong 
quantizing magnetic field, in the following form: 



p{x,x') = '^p{x,x',py,p;,,iy)exp[i{{t-t')E_ - {y - y')py - {z - z')p;,}] 



(9) 



where the sum is over the momentum components Py, Pz and the Landau quantum number v. Since the momentum 
variables are continuous, the sum over momentum components will be replaced by the corresponding integrals. Now 
we have from eqn.(9) 



+1/2 

P{X,X' ,Py,Pz,l')= ^ V{X,Py,Pz,V)v\x' ,Py,Pz,V) 

spin=—\l2 

Then substituting the negative energy up and down spinors states, we have 

1 



p{x,x' ,py,pz,y) 



2E. 



■ [E-A - PzjzloA + m^oA - pi_^y^(yB]9{K- \pz\) 



(10) 



(11) 



(see Appendix for detail derivation) where the matrices A and B are given by 

f hll 

/.-i/^_i 

ui 

\ /.-i/^_i 



(12) 



/ lu-ili 

o_ ui_, 
^0 

\o i,x_^ 

where the primes indicate the functions of x' . Now in the evaluation of vacuum energy, we have noticed that it would 
be more convenient to define a quantity pf, similar to the chemical potential for the /th flavor in a multi-quark 
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statistical system in presence of strong quantizing magnetic field (strictly speaking we are not considering a multi- 
quark statistical system and /i/ is therefore not the quark chemical potential. However, its minimum value should be 
m and not zero, i.e., in this simplified model, just like dynamical mass m, this quantity is also treated as a parameter 
and we evaluate numerically ji / and m and then obtain the upper limit of Landau quantum number [vmax] and the 
cut of A). Then it is very easy to write 

A= (M/'-m2-2i/g/B„)'/' (14) 

Since A > 0, it is also possible to express the upper limit of v, which is the maximum value of Landau quantum 
number of the levels occupied by /th flavor, and is given by 



max 



(15) 



where [ ] indicates the nearest integer but less than the actual number. Now to obtain the energy density of the 
vacuum, we consider the NJL (chiral) Hamiltonian, given by 



H = Y,m + \Y.^{z,3) (16) 

^ 1 

= ^75(»)S(*).(p; - 2(25)E^(^^ -^^Ol/JW/JQ-) -/3W75(*)/3(j)75(j)] (17) 

i=l ijkj 

where S, 75 and j3 arc usual 4x4 matrices and 2g is the effective coupling. Here we have used the formulation 
of da Providencia et al [24] of the mean field density matrix to describe the Dirac vacuum, thereby employing the 
Thomas-Fermi semi-classical method instead of formal field theory. As we have noticed, the Physics of condensation 
energy is more transparent in this method than the formal field theoretic technique. Assuming the magnetic field B„ 
along ^-direction and is constant, we can choose the gauge = {{),Q,xBm,^)- The energy of the vacuum is then 
given by 

^v= I dxitn[{-i^f..{pi - qfA)}pp,]+e^^^ (18) 

where is given by eqn.(ll), the first term is the kinetic energy part and ei^^ indicates the interaction term, including 
the exchange interaction. To evaluate the vacuum energy, we first calculate the kinetic energy term in eqn.(18). This 
quantity is proportional to the trace defined as Tr(/9/i), can easily be evaluated by using p from eqn.(ll) and the 
single particle Hamiltonian h from eqn.(17). Now using the ortho normality relations for the Hermite polynomials at 
the time of evaluation of integral over dx and also using the anti-commutation relations for 7-matrices, we have the 
first term at zero temperature (see Appendix) 

e(°) = 2A^e ^ ^ ^ (2 - 5.0) / dp.E (19) 
/=„ ^'^ v=o Jo 

where = + 2pqfBm, Nc = 3, the number of colors, and E- = —E^. 

In the evaluation of all the traces in this paper we have used the following important relation: 

MA2..B^B2..) =Tr{A^A2..B^B2..)g^^\ (20) 



Tr(7''7V7"^i^2..-BiB2..) = Tr{A^A2..B^B2..){g''" g''^ - g^^g"^ + g^^'g''^), (21) 

Tr(product of odd 7s with A and/or B) = etc. The other interesting aspects of A and B matrices are: 

i) fci^fc2A'Tr(AiA2) = {E1E2 - ki,k2z)TT{AiA2) 

ii) A:i^fc2^Tr(SiS2) = ki±.k2±Tr{BiB2) 

iii) fci^fc2'^Tr(AiB2) = kl^k''^'Tv{BlA2) = 

iv) pi^fci''p2.fc2-TV(AiB2) ¥=0= {E^.E,'^ - Pizku)P2±-k2±Tr{AiB2) 
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These set of relations are very recently obtained by us [25]. Since 7 matrices are traceless and both A and B matrices 
are diagonal with identical blocks, it is very easy to evaluate the above traces of the product of 7-matrices miiltiplied 
with any number of A and/or B, from any side with any order. 

To evaluate the interaction term, we first consider the direct part which is proportional to Tr{(5pp^)Tr{(}pp^) and 
it is very easy to show that Tr{Pj5p) = (see Appendix). Then using the orthonormality relations for Hermite 
polynomials and the anti-commutation relations for the 7-matrices, we have the direct term 



V^i, = -4gm^[V{A,m)]^ 
(The four fermion coupling is included in V, it is ~ V{1, 2)pip2) where 



(22) 



V(A,m) = ^J2 <'JB^ E (2 - ^.0) 



27r2 



v=0 



dpz 



(pi + m.,f)2 



(23) 



where TOj/J' = (m^ + 2vqfBmY^'^ (see Appendix for derivation). 

To evaluate the exchange term, we first calculate Tr{{(3pp^){(}pp^)). Now 



/3p = ^ [E-I3A + p^A^, + mA- p±B-iy] 



(24) 



Then to obtain the trace of the product of /3pi and (ip2 it is very easy to show that only the direct product terms are 
non-zero whereas cross product terms do not contribute. Therefore 



(3pi(3p2 



1 



- [EiPA + pizAj^ +mA- pi±Bjy] 



[E2PA' + p2zA'j, + mA' - p2±B'jy] 



(25) 



Then using(the orthonormality relations for Hermite polynomials) at the time of integration over dxi and dx2, the 
above trace reduces to 



{Ppipp2)dx = 



1 



[4:EiE2 + ipizP2z + 4m2] = 



^ _^ PuP2z _^ 



El E2 El E2 



(26) 



where both Ei and E2 are negative. Then in the energy contribution, after integrating over piz and p2z, the first 
term gives 



^Y^IfBm 5^ (2 -MA 



Similarly the contribution from second term is given by 



(27) 



2n2 



/2 



v=0 



(28) 



and finally, the third term is given by 



27r2 



v=0 



A +(A2+ ^2^^)1/2 



(29) 



(see Appendix for derivation of these expressions) 

To obtain the next term in the exchange part, we evaluate the trace Tr{{(3^^pp^){(3^r^pp^)), which unlike the direct 
case, gives non-zero contribution. Using the anti-commutation relations of 7-matrices and as usual with the help of 
orthonormality relations for Hermite polynomials, we finally arrive to the following result 



, 1 1 

E1E2 E1E2 \Ei E2 



^ _^ PizP2z _^ m 



(30) 



5 



(see Appendix for derivation) The contribution to the interaction energy will again be obtained if we integrate over 
Piz and j32«(done in similar manner as have been done for direct case). Then the first term is given by 



27r2 



(31) 



i/=0 



The second term is given by 



/2 



(32) 



i/=0 



The third term is given by 



m 



27r2 



and finally the fourth and fifth terms, which are identical, given by 



A +(A2+ ^2^^)1/2 



m, 



N. 



,,(/) 



m 



^^IfBm ^(2-Mln 



f=u 



u=0 



(33) 



(34) 



v=0 



Then combining all these terms we finally obtain the vacuum energy density. Since the mass m, which is assumed 

to be same for both u and d quarks, is generated dynamically, we obtain this quantity by minimizing the total vacuum 
energy density with respect to m, i.e., by putting dey/dm = 0. Simplifying this non-linear equation, we finally get 

dCy 



= -P + 2gQR = 
am 



where 



^ = ^E«/^- E(2-M 



27r2 



f=u 

d 



v=0 

\„U) 1 



2to-'^A 



ml J (A2 +m2 )V2 



2mX 



f=u 



N, 



i/=0 



X 



m 



A 



mlj (A2+m2_/)V2 



^ = ^ E E (2 - -^.o) [A - 4mX] 



f=u 



with 



i/=0 



X = In 



A +(A2+ ^2^^)1/2 



m, 



(35) 

(36) 
(37) 
(38) 

(39) 



It is therefore obvious from eqn.(35) that the trivial solution m = is not possible in this particular situation, or 
in other wards, the gap equation given by 



m = 4g'Vm 



(40) 



can not exist. On the other hand in a non-magnetic case, or for the magnetic field strength less than the quantum 
critical value, eqn.(35) reduces to the gap equation as written above (eqn.(40)). Here V is the overall contribution of 
interaction terms. Hence it is obvious that m = 0, the trivial solution exists in this non-magnetic or the conventional 
scenario, investigated by Bhaduri et. al. [23]. The phase with m = is the Wigner phase and m 7^ is the so called 
Goldstone phase. Now eqn.(40) further gives 



4.gV = 1 



(41) 



which is nothing but the well known gap equation used in BCS theory. The gap equation therefore does not exist in 
presence of strong quantizing magnetic field if the Landau levels for u and d quarks are populated. 



6 



3. CONCLUSION 



The non-existence of trivial solution (m = 0) indicates the spontaneously broken chiral symmetry in presence of 
strong quantizing magnetic field. Therefore as soon as the Landau levels are populated for light quarks in presence of 
strong external magnetic field, the chiral symmetry gets broken, the quarks become massive and the mass m (assumed 
to be same for both u and d quarks) is generated dynamically. 

Therefore we may conclude here that the Wigner phase does not exist in the case of relativistic Landau dia-magnetic 
system. Further, if the deconfinement transition and restoration of chiral symmetry occur simultaneously, or in other 
wards, if the c;liiral symmetry remains restored within the bag, as it is generally assumed, then it puts a big question 
mark whether the idea of bag model is applicable at all in presence of strong quantizing magnetic field. Questions 
may also arise, that if Wigner phase still exists inside the bag, then whether the external quantizing magnetic field 
c;an penetrate the bag bomidary, if not, what is the underlying physics which prevents the external magnetic field 
from entering the periphery of the bag. 

Now to illustrate the variation of dynamical quark mass with magnetic field, we consider the relation 

ml = -'^< ^i' > (42) 

where rn^ is the pion mass, mo is the quark current mass and /,r is the pion decay constant. Using the spinor solutions 
given by eqns.(6) and (7) we get 

2 2mom -^c . ^ -, 

= —f2- (2 - 5,o)ln 

We have now solved eqns.(35) and (43) numerically to obtain A and m for various values of magnetic field strength. 
In the actual numerical work we have solved self-consistently for the dynamical mass to and the parameter [if using 
eqn.(15) for z^^^x ^"^^ then from eqn.(14) we get the infra red momentum cut ofi^ A. In our calculation we have 
always used /u/ instead of A which allows us to obtain vm.wJ . So we can not compare our result with those obtained 
with zero chemical potential, since in our calculation it is just a parameter, it has no resemblances with chemical 
potential in a finite density quark (m, d) matter system. We were forced to make this trick to obtain m. A and also 
^max self-consistently from the numerical solutions of eqns.(35) and (43). In doing numerical calculations, we have 
considered the following sets of numerical values for the parameters. The current quark mass toq = lOMeV, pion mass 
TOtt = 140MeV, pion decay constant = 93MeV, coupling constant g = lOGeV^^ and electron mass m,. = 0.5MeV. 
In fig.(l) we shown the variation of dynamically generated quark mass with the strength of magnetic field. As it is 
evident that the dynamical quark mass never goes to zero and diverges beyond Bm ~ lO^^G. 

1 I 1 1 iiii| — I I 1 1 iiii| — I I 1 1 iiii| — I I 1 1 iiii| — I I 1 1 nil 
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FIG. 1. The variation of dynamically generated quark mass with the strength of magnetic field (expressed in 
BW(e) ^ 4 4 X iqISg.) 



APPENDIX A 



Eveduation Of Density Matrix: 

To obtain the density matrix (eqn.(ll)), we use eqns.(6) and (7), then 



1 



2E{E - m) 



( Pzlu \ 

-i{2vqfBmf/^I,^i 



\ 



{E - m)h 




(pX i{2uqfBm)'/^ll_, {E-m)ll o) 



and 



2E{E - to) 

-ip,{2vqfBmY'^h_X 
p,{E-m)hIl 



ip,{2vqfB^f/^lX_^ 
2vqfBmIv-iIl-i 
i{2vqfBmY'\E-m)Ul_^ 




p,{E-m)I,Il Q\ 

-i{2uqfB^f/^(E-m)h.iIl 

{E - mfhil 

0/ 



2E{E - to) 



/ i{2vqfBmY/^I,_i \ 

-Pzh-\ 



{E - to)/,_i 



(-i{2vqfBmY'^I, -pX-1 {E-m)ll_^) 



2E{E - to) 

/ (2vqfBm)hIl 
+ip,{2vqfBmY'H,_Jl 

\-i{2vqsBmf/^{E-m)h-iIl -p^{E - m)h-ill_^ 



ip,{2vqjBrnY'H,Jl_^ i{2vqfBmf/^{E -m)hll_^\ 
p%-ill^, -{E - m)pj,-ill_, 


{E - m)H.-iIl_, J 



Adding 



^T^Tt _|_ ylyl^ 



1 

2E 



( {E + m)Ul pJX i{2vqfB^)^l^hIl_^\ 

(£; + to)/,_i/^_i -i{2vqfB^Y/^h_Jl -pJ^_Jl_, 

PzUl i{2vqfB^f/^UIl_^ {E-m)Ul 

\-i{2vqjBrry'^{E-m)h^iIl -pj„^ill_^ {E - m)I,-iIl_J 

which may easily be simplified after a little algebra and reduces to 

1 



p{x,x' ,py,pz,u) 



2E- 



[E_A - PzlzloA + m.'joA - p_L7j,7oB]^(A- | Pz |) 



where the matrices A and B are given by eqns.(12) and (13). 
Evaluation of Kinetic Term: 
Free Hamiltonian 

h = 7^S.p + Pm 

substituting E from eqn.(4) and using the properties of 7 -matrices we get 



h = d.p + (3m (A8) 
the usual one. With the gauge choice = (0, 0, xBm, 0), we have 

h = 75S.(p - QfA) + 13m (A9) 

by the same technique as above, we have 

h = d.{p- qfl) + (3m (AlO) 

Now to evaluate Tr (/ip), where 

h = (7.(p - qfl) + m)7o (All) 

we take Ay = qfxBm, whereas all other components are zero. Then we can express the above Hamiltonian in the 
following form 

h = (jxPx + lyiPv - qfxBm) + IzPz + m)jo (A12) 

substituting 

(a/S™)^/^(^-x)=-C (A13) 

we have 

h = {ixPx - {(IfBruf'^^C + ^zPz + m)^o (A14) 

In evaluating the trace Tr(/i/7) wc integrate over x-coordinate, use orthonormality relations for Hi,{C,) and finally using 
the anti-commutation relations for 7-matrices, we have the first term of the product oc (from the conclusion drawn 
just after eqn.(21)), second term oc 4p^, third term oc 4m^ and the fourth term a Ap\. Finally adding, we have 
Tr(pft,) cx 2E . Since E <Q, this is actually oc —2E, where E = {p1+p']_ + m^)^/^ Then the kinetic energy density is 
given by 



4°^=2A^cE^ E(2-M/ dp^E 

f=u u=0 -"^ 



\,Af> 

N, 



f=u i/=0 



+ ml J In 



A+(A2+m2^)V2 , 



(A15) 



m^j 

where m^j = + 2vqfBm and e„ = 2/3e, = l/3e and Nc = 3 
Interaction Energy (Direct): 

It is very easy to show that Tr(/J75p) = 0: 

/375P = Ixlvlz ^ [EA + p^Aj, + mA- p±B'yy] (A16) 

Then it is very easy to check that the above result follows from here. 

To calculate the direct term of interaction we have to evaluate the Tr(/3p) integrated over a;-coordinate, i.e., 

Ti{(3p)dx = — / [E(^A + p,A-f, +mA- p^B^y]dx (A17) 

-cc J~oo 

since first, second and fourth terms contain odd (single) number 7 with A the contribution of those terms are zero. 
Only third term contributes to the integral and is given by 

— y ^ Tv{A)dx = ^J ^ 2[Ir.u'{x)+ll_,{x)]dx (A18) 

using orthonormality relation for Hermite polynomials the above integral reduces to 2m/ E. Since there is an identical 
expression with integral over x', finally we get 

ydir = -45m2[V(A,m)]2 (A19) 
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Momentum Integral: 
First Term: 



27r2 



/=u j/=0 



Second Term: 



27r2 



.,(/) 



27r2 



I maxJ 



P2zdp: 



'2z 



Third Term: 



2n^ 



N. 



A +(A2+ ^2^^)1/2 



m: 



Interaction Energy (Exchange Terms): 

The product 



Now 



{015 Pi) {015 P2) = lil2l3w^[EiAi +puAi73 + mAi -^1^^172] x 

717273 :r4r [-^2^2 +P2zA2j3 + mA2 - P2±B2j2] 

Zhi2 



717273 7^Tr[-E'i^i +PuAi73 + ^^1 -pi±Si72] 



2^1 



[-^1717273^1 +Pu7i72^i +7i7273TO>li - Pi±7i7273-Bi72] 



Hence 



(/375Pl)(/?75P2) 



1 



-[i^i 71 7273^1 +Piz7i72^i +7l7273»T^^l - Pi±7i7273-Bi72] 



AE\E2 

[-£^2717273^2 + P2z7l72^2 + 7l7273?Tl^2 - P2±7l 7273-0272] 
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Now we use the results 



717273717273 = -717273717372 = 71727172 = -1, 71727172 = -1 



(A26) 



to obtain the Tr[(/375/9i)(/375P2)], which is given by 



Tr[(/375Pi)(/375P2)] = 



1 



[-AE1E2 - '^PizP2z - 4m^ - 4mi?2 - 4m£'i] 
PlzP2z , ™^ , 1 , 1 A" 



AE1E2 



= - 1 + 



El E2 El E2 \ El E2 




(A27) 



The integration over piz and p2z are done in the same manner as are did for Tr[(/3pi)(/3p2)]- 
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